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摘要 :我 国 探 月 工程 嫦娥 七 号 将 构建 由 中 继 星 上 4. lm 口径 抛物 面 天 线 及 地 基 射 电 望 远 镜 
组 成 的 首 个 月 地 空间 其 长 基线 干涉 测量 (Very Long Baseline Interferometry, VLBI) 试验 
系统 。 不 同 于 地 基 甚 长 基线 干涉 测量 望远镜 ， 受 中 继 星 轨道 扰动 和 设备 时 延 影响 ， 射 电源 预 
报时 延 模型 和 实际 时 延 有 时 可 能 会 有 较 大 差异 ， 导致 无 法 引导 VLBI 相 关 处 理 机 正常 工作 , 需 
要 通过 射电 源 条 纹 搜索 方法 寻找 实际 时 延 模型 才 角 EE 解决 这 一 问题 ,由 于 4. Im 空间 天 线 口径 较 
小 ， 接 收 信和 号 能 力 弱 。 为 提高 信 噪 比 ， 需 要 长 时 间 积 分 ， 而 线性 时 延 模型 无 法 保证 长 时 间 的 
时 延 精 度 。 针 对 这 一 问题 ， 提 出 了 一 种 用 于 空间 甚 长 基线 干涉 测量 的 射电 源 条 纹 搜索 算法 ， 
首先 通过 建立 时 延 和 时 延 率 的 二 维 搜索 网 格 ; 使 用 每 个 网 格 点 对 应 的 先 验 时 延 和 时 延 率 构建 
线性 时 延 模型 ; 再 使 用 该 线性 时 延 模 型 对 原始 数据 分 时 段 相 关 处 理 ; 然后 采用 二 维 傅 里 叶 变 
换 搜 索 残 余 时 延 和 时 延 率 , 确定 最 大 的 相关 幅度 值 对 应 的 网 格 点 ; 最 后 使 用 该 网 格 点 的 先 验 
时 延 和 时 延 率 以 及 残余 时 延 和 时 延 率 重 构 二 次 项 时 延 模型 。 通 过 对 RadioAstron 空 间 其 长 基 
线 干涉 测量 的 地 空 基 线 实际 数据 处 理 ， 验 证 了 该 射电 源 条 纹 搜 索 算法 的 正确 性 。 
关键 词 : ”空间 VLBI; 条 纹 搜索 ， 时 延 模型 ， 相 关 处 理 ;， 射 电源 
中 图 分 类 号 : TN216 ”文献 标识 码 : A 文章 编号 : 1672-7673 (2022) 06 


其 长 基线 干涉 测量 是 一 种 具有 高 分 辩 率 、 高 精度 特点 的 射电 干涉 技术 "。 它 采用 原子 钟 
a A ARRIGO, HES MERRIE, ES 
据 并 发 送 到 异地 的 数据 处 理 中 心 , 经 处 理 得 到 观测 结果 。 在 VLBI 观 测 中 ,两 台 射 电 望远镜 间 
的 距离 称 为 基线 ， 它 决定 甚 长 基线 干涉 测量 的 最 大 分 辩 率 品 ， 若 有 空间 望远镜 与 地 面 望远镜 
构成 空地 甚 长 基线 干涉 测量 系统 ， 基 线 长 度 超过 地 球 直径 ， 可 实现 更 高 的 分 辨 率 。 我 国 探 月 
工程 四 期 任务 将 利用 嫦 娥 七 号 中 继 星 上 的 4. Im 口径 抛物 面 天 线 , 进行 月 地 空间 其 长 基线 干涉 
测量 试验 。 不 同 于 常规 地 基 甚 长 基线 干涉 测量 所 ， 受 中 继 星 轨道 扰动 和 设备 时 延 的 影响 ， 射 
电源 预报 时 延 模型 可 能 和 实际 时 延 有 较 大 差异 , 无 法 引导 相关 处 理 机 正常 工作 。 通 过 射电 源 
2 条 纹 搜索 ,可 确定 实际 相关 处 理 机 所 需 的 高 精度 时 延 模型 ,用 于 空间 甚 长 基线 干涉 测量 空地 
基线 相关 处 理 。 
目前 ， 在 探 月 工程 中 ， 探 测 器 信号 相关 处 理 时 ， 采 用 差分 单 向 测 距 (Differential 
One-way Ranging，DOR) 信 号 的 点 频 信号 特征 来 设计 条 纹 搜 索 算法 “"。 而 射电 源 信号 是 白 
噪声 信和 号, 无 法 利用 现 有 探测 器 信和 号 搜索 算法 。 目 前 针对 射电 源 信号 的 条 纹 搜索 有 eFFT 条 纹 
搜索 算法 、CAF-W 条 纹 搜索 算法 ”和 多 重 网 格 条 纹 搜 索 算法 ”。eFFT 条 纹 搜 索 算法 和 CAF-W 条 
纹 搜 索 算 法 可 实现 残余 时 延 和 时 延 率 的 计算 , 但 主要 用 于 强 源 的 时 延 模型 搜索 。 多 重 网 格 条 
纹 搜 索 算 法 的 基本 思想 是 建立 多 层次 、 不同 尺寸 的 网 格 ,通过 这 些 网 格 , 在 不 同 范围 内 搜索 
需要 的 目标 值 。 多 重 网 格 条 纹 搜索 算法 使 用 线性 时 延 模型 进行 相关 处 理 ， 在 时 延 -时 延 率 平 
面 中 搜索 条 纹 , 适用 于 信号 较 强 , 短 时 间 积 分 就 能 得 到 条 纹 的 应 用 场景 。 由 于 空间 小 口径 天 
线 接收 信号 较 弱 , 获得 干涉 条 纹 需要 进行 较 长 积分 时 间 。 使 用 通常 的 线性 时 延 模型 无 法 保证 
较 长 时 间 范 围 内 的 相关 处 理 机 模型 时 延 精 度 ， 需 设法 构造 精度 更 高 的 二 次 项 时 延 模型 。 如 果 
常规 多 重 网 格 中 增加 时 延 二 次 项 , 会 使 搜索 空间 从 二 维 扩展 到 三 维 , 计算 复杂 度 增加 一 个 
量 级 。 
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网 格 条 纹 搜 索 算 法 , 提出 了 一 种 用 于 空间 甚 长 基线 干涉 测量 的 射电 源 条 纹 
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搜索 算法 : 首先 建立 关于 时 延 和 时 延 率 的 二 维 搜索 网 格 , 使 用 每 个 网 格 点 对 应 的 先 验 时 延 和 
时 延 率 信 息 构 建 线性 时 延 模型 ,对 原始 数据 分 时 段 相 关 处 理 ; 然后 由 二 维 传 里 叶 变换 求解 残 
余 时 延 和 时 延 率 后 进行 补偿 , 得 到 每 个 时 段 的 时 延 和 时 延 率 , 计算 每 个 网 格 点 对 应 的 相关 幅 
度 值 ; 最 后 在 最 大 相关 幅度 值 对 应 的 网 格 点 上 重 构 二 次 项 时 延 模 型 。 采 用 该 算法 ， 对 俄罗斯 
RadioAstron 空 间 其 长 基线 干涉 测量 的 实际 观测 数据 进行 相关 处 理 时 取得 了 较 好 的 效果 。 


1 射电 源 条 纹 搜 索 算法 


目前 中 国 甚 长 基线 干涉 测量 网 软件 相关 机 使 用 五 次 项 时 延 模型 进行 相关 处 理 。 在 条 纹 
搜索 的 过 程 中 , 为 减少 计算 复杂 度 , 通常 使 用 由 时 延 常 数 项 和 时 延 一 次 项 系数 构成 的 线性 时 
延 模型 。 但 随 着 时 间 范 围 的 扩大 ， 线 性 时 延 模型 的 误差 与 五 次 多 项 式 模型 的 误差 也 会 增加 。 
综合 考虑 到 条 纹 搜 索 计 算 量 与 模型 精度 因素 , 我 们 设法 构建 由 常数 项 、 一 次 项 系数 和 二 次 项 
系数 构成 的 二 次 项 时 延 模型 。 
1. 1 相关 处 理 算法 

假设 两 台 站 接收 的 目标 源 射频 信号 分 别 为 s1(t) 和 s2(t)， 同 一 波 前 信号 到 达 两 台 站 时 延 
差 为 T， 则 


— 


S(t) = si1(t + T). 


日 射频 信号 得 到 的 基 频 信和 号 


jami 


Xx1(t) = s (t)e Tht, 
X(t) = ss(t)e Pt = s (bere hot, 
其 中 ， 太 为 通道 的 天 空 频率 。 
相关 处 理 以 台 站 1 为 参考 ， 使 用 预报 模型 re 对 台 站 2 进行 时 延 补 偿 ， 即 
Xs(t—T) = s(t —T, Je P™E-Te) = s (t +T — Te)e holtre), 
经 过 条 纹 旋转 后 


X2a(t 一 To)e-72r1ore = s(t —T, Je TPt = s (t +T- to)e Tht, 


Ht 


对 台 站 1 和 人 台 站 2 的 信号 进行 傅 里 叶 变换 


FFT 
x(t) = Sif + fo) 
xo(t — Tee rhore E S, (f + fo)e PMU tote = S (f + foje TOHE, 
Sif + 万) 为 台 站 1 信号 经 过 频率 转换 后 的 功率 谱 , KIRA Bk) OM 
Pralf) = [Sif + fo)[ ee 


由 《7) AA, MRR ET, 一 t 达 到 最 小 值 时 ， 两 全 站 的 互相 关 功 率 谱 达 到 最 大 峰值 。 
1.2 残余 时 延 搜索 算法 


a) 
(2) 


(3) 


(4) 


(5) 


(6) 
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Correlation amplitude 


Residual delay rate 2 -1 f Residual delay 


图 1 2-D 傅 里 叶 变 换 残余 时 延 搜 索 图 
Fig. 1 Search function 
根据 二 维 傅 里 叶 变换 的 性 质 可 知 ， 相 关 处 理 得 到 的 互相 关 谱 二 维 传 里 时 变换 后 ， 可 以 

得 到 一 个 关于 残余 时 延 和 时 延 率 的 二 维 互 相关 函数 "， 如 图 1， 其 公式 为 


F (At, At) = If ii Ps(f, O a 


其 中 ，At 为 残余 时 延 ，Ai 为 残余 时 延 率 ;A 为 通道 带宽 ， 7 为 一 组 二 维 傅 里 叶 变 换 的 数据 时 
长 ， 包 为 天 空 频率 ，F(AT, 人 1) 为 不 同 残 余 时 延 和 时 延 率 下 的 相关 幅度 值 ，Pi2(f,t) 为 相关 处 
理 后 的 互相 关 功 率 谱 。 

将 《〈8) 式 离 散 化 ， 即 


2BT = = -i -i 
| 


其 中 ，F(ay,b,) 为 经 过 二 维 傅 里 叶 变 换 后 得 到 的 相关 幅度 值 ”;， a 和 b, 分 别 为 离散 化 后 的 残 

余 时 延 和 时 延 率 ; 人 为 频 域 的 傅 里 叶 变 换 点 数 ，/ 为 时 域 的 傅 里 叶 变 换 点 数 。 

1.3 条纹 搜 索 与 二 次 项 时 延 模型 重 构 流程 
根据 以 上 原理 , 对 于 空间 其 长 基线 干涉 测量 射电 源 条 纹 搜索 与 二 次 项 时 延 模型 重 构 流程 

如 图 2。 


(8) 


(9) 
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(Start) 
Constructing 2D search grid 


wd.。  -------------------------------------- 寺 -------------------------------------- rr------------|------------ 


:| 2D grid points(1,1) |: : 2D grid points(i,j) : :| 2D grid points(m,n) : 
: : : Correlation processing : : : 
Divide two-dimensional FFT groups 


FFT group | : : 


Two dimensional 


: ; Two dimensional 
: FFT group k 


FFT group 1 


Search for residual : 
:| delay and delay delay and delay delay and delay Į: 
hk rate rate rate : 


:ff Calculate the |] : : if Calculate the |; 
:| average correlation |: : : :| average correlation |: 
: amplitude : : Calculate the average correlation amplitude : H amplitude : 


:| Search for residual Search for residual 


Find the grid point 
corresponding to the 
maximum correlation 

amplitude 


term coefficient 
图 2 射电 源 条 纹 搜索 搜索 与 二 次 项 时 延 模型 重 构 流 程 图 


Fig.2 Flow chart of radio source fringe search and quadratic term delay model reconstruction 


由 二 维 侍 里 叶 变换 的 性 质 可 知 ， 相 关 后 二 维 伟 里 叶 变 换 的 残余 时 延 和 时 延 率 的 范围 分 
aal- gh [= 坊 , 直 ， 因 此 ， 时 延 和 时 延 率 两 个 维度 搜索 步 长 应 在 此 范围 内 。 根 据 


4B’ 4B. 2foT’ 2foT 
对 实测 数据 的 处 理 经 验 ， 设 置 时 延 步 长 为 
dstep = x35? (0) 
时 延 率 步 长 为 
dr r ay 


条 纹 搜 索 的 步 长 与 时 延 模型 精度 和 计算 量 密切 相关 ， 如 果 搜 索 步 长 设置 过 大 ， 可 能 得 
到 的 模型 精度 太 低 无 法 用 于 后 续 数据 处 理 ; 如 果 步 长 设置 太 小 , 由 于 目前 的 算法 采用 全 网 格 
搜索 策略 , 依次 计算 每 个 网 格 点 数据 会 极 大 影响 整体 的 计算 效率 。 注意 到 步 长 设置 较 小 时 可 
能 会 有 多 组 网 格 点 搜索 到 满足 精度 要 求 的 模型 , 后 续 会 研究 基于 二 维 随机 游 走 的 网 格 搜索 策 
略 ， 找 到 合适 的 时 延 模型 就 停止 计算 ， 这 样 可 以 有 效 地 提高 效率 。 

确定 步 长 后 ， 根 据 搜 索 范 围 和 先 验 时 延 模型 划分 网 格 点 。 每 个 网 格 点 对 应 不 同时 延 和 
时 延 率 ， 代 表 经 过 不 同时 延 补偿 的 数据 处 理 。 如 图 3， 以 先 验 时 延 模型 中 的 时 延 和 时 延 率 为 
中 心 点 ， 等 步 长 计算 每 个 网 格 点 对 应 的 时 延 和 时 延 率 值 。 

由 于 弱 信 号 使 用 线性 时 延 模型 进行 长 时 间 积 分 无 法 保证 时 延 精度 ， 获 得 干涉 条 纹 ， 所 
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以 将 原始 数据 划 为 7 组 较 短 的 数据 段 进行 处 理 ， 每 一 组 为 一 个 二 维 传 里 叶 变 换 组 ， 时 间 长 度 
为 7。 在 每 个 网 格 点 上 构造 线性 时 延 对 原始 数据 补偿 后 进行 相关 人 处理， 相关 处 理 的 积分 时 间 
为 T/M， 一 个 二 维 傅 里 叶 变换 组 的 数据 大 小 为 M x N。 得 到 两 台 站 的 互相 关 功 率 谱 后 ， 利 用 
CO) 式 对 每 一 组 的 相关 谱 进 行 二 维 傅 里 叶 变 换 ， 得 到 每 个 网 格 点 上 的 7 组 残余 时 延 、 时 延 率 
和 相关 幅度 值 , 利用 最 大 相关 幅度 的 平均 值 确定 最 优 解 对 应 的 网 格 点 , 在 该 网 格 点 上 进行 时 
延 和 时 延 率 的 拟 合 ， 得 到 最 优 的 二 次 项 时 延 系数 。 


1 


一 -一 一 一 A grid point 


图 3 二 维 网 格 示意 图 


Fig.3 Schematic diagram of two-dimensional grids 

在 拟 合 时 ， 将 每 一 组 的 残余 时 延 Cf RUIN ERT, EG AED EE SHR I SE A HEN 
延 qf 和 时 延 率 bf 相 加 ， 可 以 得 到 该 组 的 时 延 Q; 和 时 延 率 b;:， 即 

a; = afta, 
b; = bf + bY. 

其 长 基线 干涉 测量 观测 具有 时 间 连 续 性 ， 时 延 模 型 也 应 该 是 连续 的 。 但 在 求解 残余 时 
延 的 过 程 中 ,数据 分 段 导 致 每 一 组 之 间 的 时 延 具 有 间断 点 ,在 拟 合 二 次 项 时 延 系数 时 ,需要 
根据 连续 的 约束 条 件 ， 即 保持 第 1 组 的 时 延 wo 不 变 ， 修 正 其 余 组 的 时 延 wi ， 使 之 与 前 一 组 末 
的 时 延 相等 ， 即 


ag = ao 
i = ai1 + biT’ 
Ep, aP NBERE; 7 为 二 维 傅 里 叶 变 换 组 的 个 数 。 
经 过 上 述 时 延 改 正 后 ， 再 进行 二 次 项 时 延 系 数 的 拟 合 。 设 未 知 的 二 次 项 时 延 系数 为 
[xo xi Xx2]， 求 解 时 延 的 多 项 式 为 


0O<i<l ， 


a” = xo + xiti + Xt, 
其 中 ，ti = i x 7 为 每 个 二 维 傅 里 叶 变 换 组 开始 的 时 间 。 
对 《15) 式 右 端 求 导 ， 得 到 时 延 率 的 多 项 式 为 
bi = Xi + 2x3ti. 
联 立 (15) RA G6 式 建立 超 定 方程 组 ， 再 利用 正 交 三 角 分 解法 求解 ”。 


Ht 


2 算法 验证 


RadioAstron 项 目 是 苏联 在 20 世 纪 80 年 代 中 期 提出 的 一 项 空间 其 长 基线 干涉 测量 
(SVLBI) 计划 。2011 年 , 俄罗斯 发 射 了 10m 口 径 的 RadioAstron 空 间 射 电 望 远 镜 (Spektr-R, 
Ra) 卫星 。Ra 与 地 面 多 个 大 型 射电 望远镜 构成 了 空间 其 长 基线 干涉 测量 系统 。 其 接收 频段 为 
0.3GHz, 1.6GHz, 5. 0GHz 和 22GHz， 轨 道 高 度 近 地 点 600km， 远 地 点 330000km， 通 过 联合 地 面 
望远镜 进行 干涉 测量 。 

为 验证 空间 其 长 基线 干涉 测量 射电 源 条 纹 搜索 算法 ， 我 们 使 用 2014 年 空间 射电 望远镜 
Ra 与 美国 阿 雷 西 博 射电 望远镜 (Arecibo Radio Telescope, Ar) 联合 观测 射电 源 0823+033 
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的 数据 进行 干涉 处 理 。 观 测 起 始 时 间 为 2014 年 4 月 11 日 23 时 07 分 00 秒 ， 结 束 时 间 为 2014 年 4 
月 11 日 23 时 07 分 59 秒 。 其 中 Ar 站 为 2 比特 量化 ，4 比 特 户 出 的 数据 ， 空 间 站 Ra 为 1 比特 量化 ，2 
比特 扇 出 的 数据 。 天 空 频 率 4836MHz， 带 宽 16MHz， 积 分 时 间 为 8.192 x 10-3s， 相 关 处 理 后 
的 数据 规格 为 60 x 120 x 8192, 利用 每 120 x 8192 个 点 的 数据 为 一 组 进行 二 维 傅 里 叶 变 换 。 
由 (10) 式 和 CLL) 式 确 定时 延 和 时 延 率 方向 的 搜索 的 步 长 为 


dstep = => = 1.02400 x 10-4s， 


—" = 1.26210 x 10-8s/s. 
2foT 


以 先 验 时 延 模型 作为 搜索 的 中 心 点 ， 时 延 和 时 延 率 搜 索 的 网 格 大 小 为 20 x 200， 其 中 时 

延 的 搜索 范围 为 6.33600 x 10-3s~8.28160 x 10-3s， 时 延 率 的 搜索 范围 为 1.62700 x 
-6s/s~4.13858 x 10-6s/s。 

算法 验证 时 ， 首 先 采 用 常规 多 重 网 格 条 纹 搜索 算法 ， 仪 使 用 线性 时 延 模 型 ， 对 每 个 网 
格 点 相关 处 理 , 得 到 峰值 点 的 干涉 条 纹 图 。 再 使 用 本 文 研 究 的 射电 源 条 纹 搜 索 算 法 构建 二 次 
项 时 延 模型 ， 最 后 对 两 者 的 结果 比较 验证 。 

图 4(a) 为 直接 利用 线性 时 延 模型 进行 条 纹 搜 索 的 相关 幅度 图 ， 其 峰值 点 对 应 的 网 格 位 
置 为 (15，158) ， 相 关 幅 度 值 为 2.86900 x 10-4， 图 4 人 b) 为 峰值 点 对 应 的 条 纹 相 位 。 图 中 
没有 明显 的 干涉 条 纹 。 
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图 4 使 用 条 纹 搜索 和 线性 时 延 无 法 得 到 正确 条 纹 。 (a〉 相 关 幅 度 图 ; O) 干涉 条 纹 图 


Fig. 4 Fringe search using linear delay can not produce the fringe. (a) Correlation amplitude; (b) interference fringe 


图 5 为 空间 其 长 基线 干涉 测量 射电 源 条 纹 搜 索 算法 得 到 的 不 同 网 格 点 对 应 的 相关 幅度 
值 ， 峰 值 点 对 应 的 网 格 位 置 为 (11，101〉， 相 关 幅 度 值 为 1.96700 x 10 飞 ， 经 过 拟 合 后 得 
到 二 次 项 时 延 系数 [7.36875 x 1073s, 2.88784 x 10-6s/s， 一 4.29989 x 10-11s/s?|。 使 用 
该 二 次 项 模型 做 相关 处 理 ， 可 得 到 图 6 (a) 的 干涉 条 纹 。 图 6 (b) 为 解 卷 绕 之 后 的 条 纹 结果 。 可 
见 ， 本 文 提出 的 算法 在 处 理 RadioAstron 数 据 时 要 优 于 多 重 网 格 条 纹 搜 索 算 法 。 
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图 5 二 次 项 时 延 补偿 的 相关 幅度 


Fig.5 Correlation amplitude using quadratic delay 
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图 6 采用 二 次 项 时 延 补偿 后 获得 清晰 的 干涉 条 纹 。 
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Frequency/Hz «108 


(b) 


(a) 解 卷 绕 之 前 的 干涉 条 纹 ;， Cb) 解 卷 绕 之 后 的 干涉 条 纹 


Fig. 6 Clear interference fringes are obtained by quadratic time delay compensation. (a)Phase wrapping; (b) phase 


unwrapping 


该 算法 原理 样机 在 Linux 集 群 平台 ， 采 用 计算 节点 间 MPI+ 节 点 内 POSIX 多 线程 的 两 级 并 
行 实现 ， 计 算 节 点 的 中 央 处 理 器 型 号 为 英特尔 的 Xeon E5-2620， 主 频 为 2GHz， 每 个 节点 的 中 
央 处 理 器 个 数 为 12， 运 行 耗 时 如 表 1。 后 续 通 过 优化 ， 运 行 速度 可 进一步 提高 。 


表 1 空间 VLBI 射电 源 条 纹 搜索 算法 


Table 1 Time consuming of the radio source fringe search algorithm for space VLBI 


运行 耗 时 


Number of nodes*Number of threads 


Grid points (delay*delay rate) 


Time consuming 


10*10 


20%200 


About 45 minutes 


3 结论 


在 将 来 探 月 四 期 任务 中 ， 地 月 空间 甚 长 基线 干涉 测量 试验 系统 受 中 继 星 轨道 扰动 和 设 
备 时 延 影 响 , 可 能 导致 射电 源 预 报时 延 模 型 和 实际 时 延 有 较 大 差异 。 本 文 在 常规 多 重 网 格 条 
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纹 搜 索 基础 上 , 根据 空间 其 长 基线 干涉 测量 信号 弱 的 特点 , 提出 了 一 种 适用 于 空间 甚 长 基线 
干涉 测量 的 射电 源 条 纹 搜索 算法 , 该 算法 先 建立 关于 时 延 和 时 延 率 的 二 维 搜索 网 格 , 使 用 线 
性 时 延 相关 处 理 ,然后 由 二 维 传 里 叶 变 换 搜 索 残 余 时 延 和 时 延 率 , 最 后 重 构 二 次 项 时 延 模型 。 
本 文 使 用 RadioAstron 数 据 验 证 算法 ， 结 果 表 明 在 使 用 常规 线性 时 延 模 型 搜索 无 法 实现 空地 


基线 条 纹 的 情况 下 , 新 算法 可 以 得 到 清晰 的 空间 其 长 基线 干涉 测量 射电 源 干涉 条 纹 。 经 进 一 
步 测 试 和 完善 ， 该 算法 可 应 用 于 我 国 后 续 空 间 其 长 基线 干涉 测量 。 
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Abstract: The fourth phase of China's lunar exploration project will build the first lunar 
terrestrial Space VLBI test system composed of 4.1m antenna on Chang'e 7 relay satellite and 
ground-based radio telescope. Different from the ground-based VLBI system, due to the influence 
of relay satellite orbit disturbance and equipment delay, the predicted delay model of radio source 
may be quite different from the actual delay, which can not guide the relevant processors as the 
core equipment of VLBI data processing to work normally. In this paper, a radio source fringe 
search algorithm for Space VLBI is proposed. Aiming at the characteristics of weak signal 
receiving ability of space small aperture antenna and the need to use quadratic delay model to 
ensure long-time integration accuracy, firstly, a two-dimensional search grid about delay and 
delay rate is established, Using the prior delay and delay rate information of each grid point, the 
linear delay is constructed to process the original data in different time periods, and then the 
residual delay and delay rate are searched by two-dimensional Fourier transform to determine the 
grid point corresponding to the maximum correlation amplitude value. Finally, the quadratic delay 
model is reconstructed by using the delay corresponding to the grid point. The algorithm can be 
applied to the radio source fringe search of the ground air baseline of the RadioAstron project. 


Key words:Space VLBI; Fringe search; Delay model; Correlation processing; Radio source 


